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ABSTRACT 


Cooling of neutron stars is calculated using an exact stellar evo- 
lution code. The full general relativistic version of the stellar 
structure equations aze solved, with the best physical input currently 
available. For neutron stars with a stiff equation of state, we find 
that the deviation from the isothermality in the interior is significemt 
and that it takes at least a few thouswd years to reach the isothermal 
state. By comparing the most recent theoretical and observational 
results, we conclude that for Cas A, SN1006, and probably Tycho, 
"f.tamdcurd" cooling is inconsistent with the results from the Einstein 
Observatory, if neutron stars are assumed to be present in these objects. 
On the other hand, the "detection" points for RCW103 and the Crab are 
consistent with these theoretical results. 

Subject headings : stars : neutron — stars ; interior — 

nebulae : supernova remnant — X-rays : sources. 
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I. INTRODUCTION 


Until recently, the problem of neutron star cooling has been mainly 
of academic interest. Indeed, the sole information on the direct 
observation of surface radiation from neutron stars has been the upper 
limit for the Crab (Toor and Seward 1977) . The exciting results from 
the Einstein Observatory are already changing the whole picture. As of 
March 1981, we have at least eight upper limits (Helfand 1981; Murray 
1981; Pye et al. 1981; Seward 1981) and nearly a dozen possible "detec- 
tions" (Harnden 1981; Helfand 1981; Tuohy and Garmire 1980) which may be 
related to neutron stars in young supernova remnants less th 2 ui a- lO^yrs 
old. Prompted by these recent observational developments, we have re- 
exaunined theories of neutron star cooling. The most important factor 
which has been neglected in most of the earlier calculations is the 
effect of the finite timescale of thermal conduction. 

In these calculations (Glen and Sutherland 1980; Tsuruta 1979; Van 
Riper and Lamb 1981) , it was assumed that tne central core is isothermal 
and that the thin outer envelope has constant photon luminosity. Under 
these assumptions, the surface temperature instantaneously follows the 
decrease in the central temperature through the escaping neutrinos . 
However, such a treatment is valid only if the timescale of thermal 
conduction is infinitely short as compared with the timescale of evolu- 
tion. In the early stages of neutron star cooling, neutrino emission 
is much faster than thermal conduction, so that the thermal behavior of 
the surface layers does not necessarily keep pace with the neutrino 
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cooling in the Interior. It has been shown that this effect is negli- 
gible for models with soft equations of state with age t > 10 yrs 

(Malone 1S74; Richardson 1980) . Ray (1981) suggested through analytic 
estimates that, for some stiff equation of state stars with a leurge 
crust, this effect may be important because such a crust has a smaller 
conductivity and contains a significant fraction of the stellar mass. 
This question can be accur^t> ly answered by solving the equations of 
stellar structure siroulte jeously as a mixed initial -boundary-value 
problem. We have made such calculations for "standard" cooling of the 
stiff model by using an exact stellar evolution code. (By "standard" 
cooling we mean, following earlier convention, cooling in the absence 
of such exotic particles as pion condensates and quarks.) We have also 
recalculated the soft model using the isothermal method, because we are 
currently interested only in stars older than 100 yrs. In both cases, 
we have made use of the best theoretical knowledge currently available. 
In this Letter , we briefly report our results with particular emphasis 
on the stiff model and compare them with observational data. A more 
detailed description and discussion will be given elsewhere. 


II. PHYSICAL INPUT 

Physical input for the exact evolution calculation is the following 
The full general relativistic equations (Thome 1977) are adopted for 
the hydrostatic equilibrium, thermal energy transport, and energy con- 
servation, These equations are integrated simultaneously as a function 
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of time and space coordinates by using a Henyoy-type stellar evolution 
code (Sugimoto 1970) . We use a stiff equation of state by Pandharipande , 
Pines, and Smith ^1976) generally called the PS model, based on tensor 
interaction:; between nucleons, for p > 4.3 x lO^'g cm“’, and the finite 
temperature e iiation of state for a mixture of electrons and ions for 
p < 4.3 X 10^'g cm“®. (p denotes the matter density.) The density at 
the boundary of the heavy ion crust and the central nucleon liquid core 
is chosen to be 2.2 x lo'‘*g cm“\ Specific heats of all fermions 
(neutrons, protons amd electrons) and heavy ions are included and the 
effects of crystalization are taken into accovint. Neutrino emissivities 
included are: modified URCA, neutron -neutron and neutron-proton brems- 

strahlung (Friman and Maxwell 1979) , electron-ion bremsstrahlung .Scyeur 
and Brown 1979) , and plasmon, pair and photo neutrino emissivisites 
(Beaudet et al. 1967). Weinberg-Salam theory with the Weinberg angle 
sin 6^^ ■ 0.23 is applied. Electron, muon and tau neutrinos are included. 

We adopt the approach of Maxwell (1979) in taking into account the 
effect of superfluid protons and neutrons on neutrino emissivities and 
specific heats. The critical temperatures, ^cr^ ' superfluid 

neutrons and protons are given by Ta)catsuka (1972) and Chao et al. (1972), 

* * 

respectively. We set the effective mass m /m * m /m •» 0.8. We use 
the Los Alamos table (Huebnar et al. 1977) for iron radiative opacities 
and work by Flowers and Itoh (1976, 1979) for conductive opacities. 

The effect of magnetic fields on opacities is taken into account by 
adopting the approach of Tsuruta (1979) . 
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III. RESULTS 


Our results axe summarized in Figures 1 and 2, for the PS model 

where: the gravitational mass is • 1.31 (the rest mass ■ 1.40No)i 

00 

the stellar radius is R « 16.1 km (the "apparent** radius R » 18.5 km), 
the central density is ■ 4.0 x lo'‘*g cm"®, the redshift correction 
factors (Thorne 1977) are e*^ « 0.75 and 0.87 at the center and at the 
surface, respectively (4> is the gravitational potential) , the mass of 
the central nucleon core is 0.8 Me, wd the thickness of the crust is 
5.4 km. 

In Figure 1, the "stamdeurd" cooling curves are shown for four 
cases: Case S with superfluid nucleons and Case N with non-superfluid 

(normal) nucleons both for the zero magnetic field, and Cases MS (super- 
fluid) and MN (non-superfluid) for the magnetic field of H * 4.4 x 10** 
gauss. In Figure 2, the evolutioneury change in the local temperature 
distribution is shown as a function of local density for Cases S and N. 

(The stages with the same number, #0-8, correspond to the same age.) 

The effects of the finite timescale of thermal conduction are clearly 
seen in these figures. The important outcome is that (i) the temperature 
distribution in the interior Deviates significantly frcmi the isotherma- 
lity for as long as a few thousand years, and that (ii) unlike the smooth 

fall which is the characteristic of the isothermal cases, the surface 

00 

temperature T^ declines as follows. 
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a) Superfluid Neutron Stars (Case S) 


Phase 1 (t ^ 10“^yrs) : The initial (#0) model with the central 

temperature ■ 1 x 10^ ®K was constructed by suppressing neutrino 

emission and choosing a long enough time step to reach the thermally 

relaxed state. Afterweurds, neutrino emission is turned on. The neutrino 

energy loss rate, dominates over the photon cooling rate, 2 

3 (e^^L )/3 m^ (i.e., e,, » le . 1) in the interior where p > n. lo^g cm“®. 
r i V ' pn ' 

On the other hand, in the surface layers. In Figure 2, the 

line for e = e , is shown for several stages. 

V ph 

The nucleon core (p > 2 x lo'‘*g cm”^) is cooled by the modifieo 
URCA neutrino emission. The neutrons in the region 10**g cm"^ p ^ 10*** 
g cm“* become *Sq superfluid. (The corresponding maximum superfluid 
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critical temperature in this region is T^“' n, 1.7 x 10* °K.) The greatly 
decreased specific heat at T << T^^^ causes the valley in this region 
(Figure 2, #1). The trough around p 't 10 g cm is due to the plasmon 


neutrino losses; the density dependence of causes the temperature 
inversion at p % 10’-10**g cm”*. As a result, the temperature peak 


appears near p 10^-10* ®g cm”’(#l). 
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The central temperature T^ decreases monotonically during this 
phase, but the temperature near the surface (p < n/ I0*g cm”*) stays 
almost constant, because a period of 10”*yrs is too short for thermal 
conduction to be effective. 

Phase I I (10”^ < t ^ 1 yr) : This phase is characterized by the 


appreciable decline cf T^ as seen in Figure 1 (##1 - 3) . This decline 
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is caused by the decrease of the temperature peak around the outer crust 

of p " 10*-10*°g cm“* through the plasmon neutrino emission. Since 

photon cooling is dominant > e^) in the surface layers where 

p < 'V lo^g cm~\ decreases simultaneously with the temperature fall 

of this boundary region. By the end of this phase, the central dips in 

Figure 2 disappear because the URCA process begins to be suppressed by 

(d) 

the proton superfluidity at T 6 x 10 K. 

<. so 

Phase III (1 yr < t 20 yrs) : The decreases in turn out to 

be very slow (##3 - 4) as compared with Phase II. Neither thermal con- 
duction nor the plasmon neutrino vunission in the outer cinist are effec- 
tive in cooling the surface layers. Thermal conduction is not efficient 
enough to inake the interior Isothermal; the temperature in the central 
nucleon core is higher them in the crust because of the suppression of 
the URCA process. 

^ 00 

Phase IV (20 yrs < t 400 yrs) ; decreases significantly again 
(##4 - 6) . The thermal conduction timescale starts becoming comparable 
with the cooling rate due mainly to the crust neutrino bremsstrahlung 
in the inner crust. Accordingly, this neutrino cooling is effectively 
transmitted to the surface. The temperature in the nucleon core, however, 
is still apprv»ciably higher than the crust temperature. 

Phase V (400 < t 4000 yrs) : The neutron star approaches a 

thermally relaxed state (##6 - 7) , while appreciable heat flow occurs 
from the hotter core to the cooler crust layers. At t 't 4000 yts, t: e 
isothermal state is almost reached. 

Phase VI (4000 yrs < t) ; The cooling is now so slow that the iso- 
thermal approximation is valid. 
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b) Non-&uperfluid Neutron Star (Case N) 


As shown in Figure 2 , the cooling in the central nucleon core is 
much faster than in Case S because the URCA process is not suppressed 
by the proton superfluidity. The inner crust layers are hotter than in 
Case S because the neutron specific heat is not suppressed by the 
superfluidity. The cooling curve in the earlier stages, which is con- 
trolled by the plasmon neutrino emission in the outer layers, is identical 

to Case S until t = 20 yrs (until the end of Phase III) . During Phase 
00 

IV, T^ in Case N is higher than in Case S because of the higher tempera- 

or 

tures In the inner crust layers. Later in Phase V, is lower than in 
Case S due to the effect of the faster URCA cooling in the core. 

Contrary to Case S, heat flows inward from the h:*tter crust to the 
cooler core during the stages approaching the thermal releu ation (Figure 
2). In Case N, the isothermal state is reached at t 1500 yrs. 

c) Magnetic Neutron Star (Cases MS and MN) 

The interior evolution for cases MS (superfluid) and MN (non- 
superfluid) are similar to the corresponding zero-field cases S and N. 

The effects of the reduced opacity near the surface under the strong 
magnetic fields begin to be significant at t = 1 yr. Afterwards the 

cooling is slower than the zero field cases, except that the final dec- 
00 

lines of T^ are steeper when photon «nission overta)ces neutrino emissions 

> || 

as the main cooling agent (t 'v, 5 x 10 yrs) . 
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d) Other Caies 


Slmiler calculations have been carried out for the PS model with 
0.4 and 0.8 and qualitatively siailaz results are obtained. 

By using the sssie physical input, %re update the earlier calcula- 
tions (Tsuruta 1980) for the soft mode?, (the BPS model; Baym et al. 1971); 
Mq <■ 1.3 and R - 8.0 kn (R >11 kn) . We used the isothermal 
approximation, which is valid for stars older than at least 100 yrs 
(see §1) . Our new results for the BPS model (see Figure 3) are essen- 
tially the same as those reported earlier (Tsuruta 1980) . 


IV. CX^ARISCXJ WITH OBSERVATION AND DISCUSSION 


In Figure 3, the preceding theoretical results are coiQ>ared with 
observational data. The solid curves refer to models which include 
superfluidity, while the dashed curves correspond to those with no 
superfluid particles. The ufqper curves represent strongly magnetized 
starn with H ■ 4,4 x lo'^qauss, while the lower curves are for the zero 
field case. The stellar mass is 1.3 M*, which is consistent with both 
current theory and observation. The crosses indicate the observed 
upper limits to the surface temperatures of various young supernova 
rfflnnants assuming the presence of neutron stars in these objects (Murray 
1981; Fabbi^ulo 1981; Seward 1981; Helfand 1981) . The circles refer to 
possible "detection" of point sources (Hamden 1981; Tuohy and Gannire 
1980) . These points are accompanied by the estimated error bars which 
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as* B^atly due to the uncertainty in the interstellar absorption. 
(Strictly speaking, the Crab "detection" should still be regarded as an 
"upper limit". However, it should be distinguished from the other upper 
limits, because unpulsed soft X*rays were indeed detected from the Crab 
pulsar while the upper limits in the other cases were determined by the 
non-detectiCT> of point sources (Hamden 1981) . ) 

We note: 

(i) Some of the observed upper limits, for Cas A, SK1006, and probably 
Tycho, are below the "standard" theoretical curries. A natural inter- 
pretation is that there are no neutron stars in these supernova remnants 
(Tsuruta 1980} , or that these star contain such exotic particles as 
charged pion condensates (Tsuruta 1979) and quarks (Iwamoto 1980 'an 
Riper and Lamb 1981) . Since Tycho is thought to be a remnant of Type I 
Supernova (SN I) , the present result is consistent with the detonation/ 
deflagration -type explosion model for SN I which disrupts the star com- 
pletely and leaves no neutron star {*ee Nomoto 1980 for a review) . 

(ii) The "detecti'Vi" points for RCW103 and the Crab are consistent with 
tltese "standard" theoretical curves. 

The "detection" point for Vela is still preliminary (Hamden 1981) . 
Very recently, possible "detection" of point sources has been reported 
for several other supernova remnants, e.g,, 3C58 and CTB80 (Becker, 
Helfand, and Szymkowiak 1981; Helfand 1981) . Further careful investiga- 
tion is required, however, before we can properly interprete the 
observational data for these sources. 

We note that our curves for the PS model are significantly higher 
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than thoM of V«n Rlp«r and Laiab (1981), while they are aonewhat lower 
than those by Glen and Sutherland (1980). For lnatar.ee, for C^ase S at 


t ” 1000 yra, our photon luminosity is 'w 6 times higher than that of 
Van Riper and Lamb, while it is a factor of 2 lowex than that of Glen 
and Sutherland. It appears that the lower boundiry of the BPS model by 
these three groups agree within the uncertainty inherent in the current 
theories of opacity and superfluic^.ity. CXir upper boundary fox the BPS 
model in higher than those of the other two groups, because we included 
the possibility that protons may very well be in a state of p-wave 
superfluidity in a high density region where p •> 't 10*^ - 5 X 10* ^g cm”* 
(Tsuruta 1979) . See Tsuruta (1981) for a more complete comparison of 
theoretical results of different groups. (We should remark here that 
in the comparison mace above, we used the liuninosity scale in Figure 3 
of the Van Riper and Lamb (1981) paper, because their BPS temperature 
scale is in error.) 

Van Riper and Lamb (1981) reached essentially the same conclusion 
a"5 ours for SN1006 and Tycho. However, for the Crab ard Cas A, the 
conclusion they present contradicts ours, as their observed points for 
the ;e sources lie above their "standard" theoretical curves. The major 
reason for this is their misplacement of these observed points (Figure 
3 of Van Riper and Lamb '’981) . 

V. CCWCLUDING REKARKS 

Cur most significant finding is that for the stiff model the effects 
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of a finite timescale of thermal conduction are quite important. For 
the stiff PS model of mass 1.3 M^, it takes 4000 yrs for Case S aoid 
1500 yrs for Case N, before the isothermal state is reached. 

Many young supernova remnants (Cas A, Kepler, Tycho, Crab, RCW103, 
and W28) have an age of a few thousand yrs or less. Therefore, for the 
stiff model, it appears that the earlier cooling calculations using an 
isothermal approximation are inadequate and we have to resort to exact 
stellar evolution calculations such as those presented here, if we eu:e 
to make relevant quantitative comparisons between theory and the greatly 
improved observational data now available from the Einstein Observatory 
and those expected from the observational progreuns planned for the neeu: 
future . 
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FIGURE C21PTI0NS 


Figure 1: 


GO 00 

The surface temperature T and photon luminosity L . (both 
observed at infinity) are shown as a fvinction of age 
(measured at infinity) for the PS model with “ 1.3 M«. 


Figure 2 : Distribution of local internal temperature as a function of 

local matter density at different evolutionary stages 
##0 - 8 for the PS model with Mq « 1.3 M* euid H « 0. Nucleon 
superfluidities are included in Case S (upper half) , while 
they eu:e absent in Case N (lower half) . 


Figure 3 : Comparison between the observational results from the Einstein 

Observatory and the theoretical curves for "steuidard” cooling 
of the PS model (upper half) and the BPS model (lower half) . 

00 O 

(left) and 1^^ (right) are shown as a function of age 
(measured at infinity) . The numbers refer to the observa- 
tional data of: (1) Cas A, (2) Kepler, (3) Tycho, (4) Crab, 

(5) SN1006, (6) RCW103, (7) RCW86, (8) W28, (9) G350.0 - 18 , 
and (10) G22.7 - 0.2. 


16 


(i-8 ■ Bja) 001 



Figure 1. 


LOG t (yr) 



LOG T (K) LOG T (K) 



LOG P (g cm"^) 


Figure 2. 


«in to 





Log t (yr) 


F iqure 3 . 






